An Orbit Fit for the Grillmair Dionatos Cold Stellar Stream 
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ABSTRACT 

We use velocity and metallicity information from SDSS and SEGUE stellar 
spectroscopy to fit an orbit to the narrow 63° stellar stream of Grillmair and 
Dionatos. The stars in the stream have a retrograde orbit with eccentricity 
e = 0.33 (perigalacticon of 14.4 kpc and apogalacticon of 28.7 kpc) and inclination 
approximately i ~ 35°. In the region of the orbit which is detected, it has a 
distance of about 7 to 11 kpc from the Sun. Assuming a standard disk plus bulge 
and logarithmic halo potential for the Milky Way stars plus dark matter, the 
stream stars are moving with a large space velocity of approximately 276 km s~^ 
at perigalacticon. Using this stream alone, we are unable to determine if the 
dark matter halo is oblate or prolate. The metallicity of the stream is [Fe/H] 
= —2.1 ± 0.1. Observed proper motions for individual stream members above 
the main sequence turnoff are consistent with the derived orbit. None of the 
known globular clusters in the Milky Way have positions, radial velocities, and 
met alii cities that are consistent with being the progenitor of the GD-1 stream. 



1. Introduction 

Data from tidal stream debris is a valuable resource for constraining Galactic struc- 
ture. In the last decade, several streams, with both globular cluster and dwarf galaxy 



progenitors, have been discovered iii the Milky Way (llbata et al 
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streams surrounding the Milky Way, a clearer picture of how our halo was built can be de- 
veloped teuUock and Johnston 2005 ). In addition, stellar streams can be used as probes of 



the Galactic grayitational potential, and thus constra ins the shape of the dark matter halo 



(llbata et al.ll2001bl : iHelmil 12004 : iFellhauer et al.ll2006l ) 



Grillmair and Dionatod (120061 ) . hereafter GD, announced the detection of a 63° cold 
stellar stream in the Galactic halo (the stream itself we refer to henceforth as GD-1, following 
GD), using stellar density counts extracted from the Sloan Digital Sky Survey (SDSS; York 
et al. 2000). This stream is extremely narrow, less than 0.25° degrees in width, which is 
less than 50 pc at their measured distances of 7.3 to 9.1 kpc from the Sun. GD therefore 
concluded that the progenitor was a globular cluster, but the progenitor remains unidentified 
and could be completely disrupted. 

In this work, we utilize newly available Sloan Extension for Galactic Understanding 
and Exploration (SEGUE; see Yanny et al. 2009) spectroscopy of stars identified along the 
stream, that are available in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), to 
constrain the orbit of the stream and search for possible progenitors. 



Data Selection 



The SDSS imaging survey jGunn et al.lll998l . l2006h has made it possible to detect faint 
Milky Way halo streams from the spatial distribution of stars because it provides very ac- 
curate multicolor p hotometry for millions of Galactic stars . Even with the we l l-calibrated 
SDSS photometry JPukugita et al.l Il996l : ISmith et all I2OO2I : iTucker et all bood : Iflogg et al. 



2OOII : iPier et al.ll2003l ). GD separated the faint structure from the backgr ound of Milky Way 



halo and disk stars only with careful filtering and smoothing techniques ( lOdenkirchen et al. 



2OOII : iRockosil I2OO2I ). To characterize the stream in more detail and to compute an orbit, 
we reanalyze the imaging data, supplemented through DR7, and then add to it all newly 
available spectroscopic observations of GD-1 stream stars available in SDSS DR7. 
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2.1. Imaging 



Using tlie analysis of GD as a guide, stars were selected from the SDSS DR7 (lAbazajian et al 



20091 ) footprint from a rough color-magnitude box restricted to blue F turnoff and upper main 
sequence stars at distances of 7 to 20 kpc from the Sun (taking into account that the turnoff 
spans more than a magnitude of absolute magnitudes): < ((? — r)o < 0.5, 18 < go < 22. 
Ma gnitudes with subscript indicate th ose which have been corrected for reddening using 



the ISchlegel. Finkbeiner. fc David (119981 ) maps as implemented in SDSS (DR7). All SDSS 
stars in the DR7 North Galactic cap footprint which are in this color-magnitude box were 
selected and plotted in an [a, 6) stellar density diagram with pixels 0.5° on a side. Examina- 
tion of this diagram by eye showed that GD-1 stream stood out weakly from the background, 
with enough contrast so that the location of several fiducial points along the stream in {a, 6) 
could be mapped. Low-order polynomials were then fit to the positions of these points. The 
lowest order best fit was of third order: 

S = -864.5161 + 13.22518a - 0.063255440^ + 0.0001009792a^ (1) 

By comparison with likely identified spectroscopic stream members we later verified that this 
polynomial matches the stream position within about 0.6° for 130° < a < 175°, and within 
about 1.0° for 120° <a< 130° and 175° <a< 220°. 

Next, 147,537 stars within ±0.5° of this polynomial fit (a conservative width is used 
to increase the signal to noise, it is not necessary to identify every star) were selected from 
135° < a < 200°. They represent the 'on-stream' data set. For a control sample ('off- 
stream'), 158,147 stars were selected within ±0.5 degrees of a curve offset by 5 degrees in 
6 from the 'on-stream' data over a similar a r ange. A Hess diagram in {{g — r)o,(7o) of the 



difference between the on and off stream data (iNewberg et al.ll2002l ) was generated and the 
results are shown in Figure 1. 

Figure 1 shows a clear faint turnoff around {{g — r)o,(7o) ~ (0.25, 19) and upper main 
sequence descending to go ~ 22.5. Other features in Figure 1 include an imperfectly sub- 
tracted thick disk residual at {g — r)o ~ 0.45 from 14 < go < 17.5 and a residual from nearby 
M dwarfs at {g — r)o ~ 1.4 and 19 < go < 23. The stars below the turnoff are concentrated 
in a relatively narrow band of color and magnitude, as expected for a stream localized in dis- 
tance from the Sun. Since the stream varies in distance from 7 to about 10 kpc for the stars 
in this figure, the actual width of the main sequence is substantially narrower than demon- 
strated here. The blue turnoff of {g — r)o ~ 0.25 suggests a lower metallicity or younger age 
than that of the spheroid, which has a turnoff of {g — r)o ~ 0.3 and [Fe/H] ~ —1.6. We 
superimpose with black dots in Figure 1 a fiducial sequence from the cluster M92. The M92 
sequence was calculated by starting with the fiducial sequence in undereddened r, {g — r) 
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from I An et al.l (120081 ) . converting it to M„ „ , (q — r)o using E(B-V) of 0.02 and a distance 



to M92 of 8.2 kpc, as tabulated by iHarrid (jl996l ). and then shifting it to the approximate 



distance to the GD-1 stream, which is approximately 9 kpc over this range of RA (distanc e 



modulus m — M = 14.76). This low-metallicity cluster, with [Fe/H] = —2.3 (jHarrislll996l ). 
fits the main sequence and the turnoff reasonably well. We do not see the giant branch or 
the horizontal branch (which should be at go ~ 15.5) in this figure, but given the faintness 
of the stream that is not unexpected. 

We now refine the rough color-magnitude box used to detect the stream in equatorial 
coordinates. The refined box is selected to allow for a stream which changes distance by 
~ 30% over 90° < a < 260°. The box is defined as the union of three selection regions 
heavily outlined in Figure 1: A) 0.15 < (5- - r)o < 0.31,17.75 < 5-0 < 19.7 (turn off) B) 
0.15 < (^ - r)o < 0.4, 19.7 <= go < 20.5 (lower turn oflF) and C) 20.5 < ^0 < 22, 17.89 < 
go — 6.52((7 — r)o < 19.52 (upper main sequence). All stars with images in the SDSS DR7 
Northern Galactic Cap region which meet these criteria (and have {u — g)o > 0.4 to exclude 
quasars), are selected and plotted in (a, 6) in the upper panel of Figure 2. The GD-1 stream 
is faintly visible, running in an arc from about {a, 6) = (220°, 58°) through (164°, 48°), then 
down to (139°, 22°) where it crosses in front of the Sagittarius stream, and then it is not 
clearly visible as it is lost in Monoceros and other Galactic halo stars near (126°, 0°). The 
lower panel of Figure 2 presents the same data in a Galactic polar projection. The GD-1 
stream is clearly visible. There are several features in the density of stars along the stream 
of unknown origin. They could be the remains of a nearly dissolved progenitor, the result 
of interactions between the stream and the potential of the Milky Way's disk and halo, or 
unassociated spheroid substructure. 

The kinematics of the stream, detailed below, reveal that it is on a retrograde orbit 
moving through the sequence of points in the order just described for the upper panel of 
Figure 2. Numerous other dwarf galaxies, streams and halo overdensities are present in 
Figure 2; these will not be discussed further here. 



2.2. Spectroscopy 



The SEGUE survey ( lYanny et al.l l2009l ) . which is one of three surveys carried out as 
part of SDSS-II, obtained spectra of approximately 240,000 Milky Way stars toward ~ 200 
sightlines that each covered seven square degrees of the sky, with an emphasis on obtaining 
spectra of fainter halo stars. While most of SEGUE's 200 observing tiles were randomly 
distributed across the SDSS imaging footprint, a few were placed on streams of known 
interest, including the GD-1 stream. 
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All S EGUE spectra Avere pr ocessed through the standard SDSS spectroscopic reduction 



pipelines (jStoughton et al.ll2002l ). from which radial velocities accurate to about 10 km s ^ for 
objects at 5^ ~ 19.5 were determined. In additi on, the stellar sp ec tra were processed through 
the SEGUE stellar parameter pipeline (SSPP) Jlee et al.l[2nn8al E]: lAllende Prieto et allbnosl ) 



in order to obtain abundance ([Fe/H]), surface gravity (log g), and other atmospheric pa- 
rameter estimates. 

We select from the SDSS-ll/SEGUE DR7 database all measured parameters of the 
12,825 spectra of stars within 1.3° of the GD-1 stream described by Eq. 1. We further 
refined the selection to exclude objects far away from the fiducial M92 curve of Figure 1 by 
requiring that they meet these color magnitude cuts: —0.3 < ((7 — r)o < 0.15, 14 < go < 17.75 
or 0.15 < - r)o < 0.4, 17.75 < ^0 < 20.5 or 0.4 < - r)o < 0.5, 15.31 < ^0 < 19.14 
or 0.5 < - r)o < 0.59, 14 < ^0 < 19.14 or 0.59 < - r)o < 0.77, 14 < ^0 < 17.21 or 
0.77 < {g — r)Q < 1.2, 14 < (7o < 15.31. The region bounded by these cuts is outlined with a 
light line in Figure 1. Most of the 4568 remaining spectra are concentrated in 3° diameter 
patches centered on SEGUE tiles, but some are part of the SDSS-I and SDSS-II Legacy 
surveys. These latter surveys targeted nearly the entire SDSS footprint spectroscopically, 
but with few and hmited signal-to-noise on stellar targets (since the SDSS Legacy survey 
primarily targets galaxy and quasar candidates). 

We show in Figure 3 the line-of-sight. Galactic standard of rest velocities, Vgsr, for each 
star in the sample, as a function of Galactic longitude. We calculate Vgsr using: Vgsr = 
RV -|- 10.1 cos b cos / + 224 cos b sin / -|- 6.7 sin b, where RV is the heliocentric radial velocity 
in km s~^ and (/, b) are the standard. Sun-centered Galactic coordinates of each star. A sine 
curve with amplitude 110 km s~^ traces an approximate locus of nearby disk stars co-rotating 
with the Sun. Spheroid stars occupy a broad range of Vgsr centered at Vgsr = 0. Seven regions 
of interest are marked along the bottom of Figure 3, indicating areas with SEGUE plates, 
where stars identified with the GD-1 stream will be selected. The positions of these seven 
regions on the sky in equatorial coordinates are also indicated with circles and numbered in 
the upper panel of Figure 2. Regions 5 and 6 were specially targeted by SEGUE with a tile 
directly on locations along the GD-1 stream. 

From examination of Figure 3, it appears that there is an excess of stars off the rotating 
disk locus at Vgsr ~ —90 km s~^ in regions 5 and 6. To confirm that these are in fact GD-1 
stream stars, we isolate the stars in regions 5 and 6 and plot their velocity histogram in 
Figure 4. 

The distribution in Figure 4 is overlayed with Gaussians for the thick disk (dispersion 
of 30 km s~^ and an offset of yU ~ 20 km s~^), and inner halo (dispersion of 100 km s~^). A 
significant peak is detected at Vasr ~ —82 km s~^ which we associate with the GD-1 stream 
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member stars. 

We next examine the metallicity distribution of stars in this velocity peak in order to 
estimate the elemental abundance of the GD-1 stream. Later in the paper we will show that 
the individual Vgsr velocities in regions 5 and 6 are 71 ± 2 and 87 ± 2 km s~^, respectively, 
so we chose a "peak" velocity range of —97 < Vgsr < —61 km s~^. 

Figure 5 shows the SSPP abundance estimates for all stars with good metallicity esti- 
mates (for a good estimate a turnoff star generally needs to be brighter than about 51 ~ 19). 
Errors on individual stars [Fe/H] are approximately 0.3 dex for spectral type F objects with 
g < 18.5. The histogram for all abundances of stars in regions 5 and 6 are plotted with a 
light line (1311 stars), those for stars in the velocity peak of Figure 4 are indicated with 
a heavy line (115 stars). The stars with velocities of the GD-1 stream are heavily biased 
towards lower metallicity stars, compared with those of the thick disk ([Fe/H] ~ —0.7), or 
inner halo ([Fe/H] ~ —1.6). 

We estimate from Figure 4 that about 30 stars in the spectroscopic dataset are from the 
GD-1 stream. To see the metallicitydistribution of the stars in the GD-1 stream, we subtract 
a scaled version of the histogram with the light line from the histogram with the heavy line. 
The scaling factor is (115-30)/(1311-30). Since the stars in the velocity selected region 
contain a smaller fraction of thick disk stars, the subtracted histogram is oversubtracted 
at high metallcities, and hkely undersubtracted at spheroid metallcitites. The mean of the 
stars in the shaded region is [Fe/H] =-1.9, but the real metaUicity of the stream is probably 
somewhat lower than this. Bins with negative counts do not appear in Figure 5. 

We now return to the sample of stars in Figure 3, and select only those of very low 
metaUicity (—2.3 < [Fe/H] < —1.65) in order to isolate stream members from the thick 
disk and halo field star populations. The low metallicity spectra with positions, colors, and 
magnitudes that make them candidates for GD-1 stream members are shown in Figure 6. 
Several velocity peaks are now clearly separated from the disk and spheroid. We now examine 
stars in each of the seven regions numbered in Figure 6 and determine their observational 
properties. 

3. Photometric Distance estimation 

We estimated the distance to the GD-1 stream at the positions of each of the seven 
regions that it overlaps using a matched filter algorithm. Wc first generated the Hess diagram 
from SDSS DR7 data from a region about 5° wide in RA and 1° wide in Dec in the vicinity of 
each plate, centered on the polynomial fit to the GD-1 stream. Then the Hess diagram of the 
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background was generated from two regions of sky with the same angular extent on the sky, 
but offset 1.5° higher and 1.5° lower in declination. The background Hess diagram (divided 
by two to correct for the difference in sky area) was subtracted from the corresponding Hess 
diagram on the GD-1 stream. The subtracted Hess diagrams are shown in Figure 7. 

Since the GD-1 stream is of quite low metallicity, we selected the globular cluster M92 
([Fe/H]~ —2.3) to compare with the observations. W e then constructed a M92 filter Hess 



diagram with the similar met hod to iGrillmairl (l2009l ). We first broaden the M92 fiducial 



sequence from I An et al.l (l2008l ) with the SDSS photometric errors. Because we do not have 
a luminosity function for M92 stars, we used the luminosity function of M13, estimated from 
SDSS survey counts vs. magnitude for stars away from the core of M13, to create the Hess 
diagram. As before, we assume the distance to M92 is 8.2 kpc. 

The M92 filter was shifted from -0.5 to 1.5 magnitudes in r in steps of 0.05 mag. For 
each shift, we cross-correlated the M92 filter with the subtracted Hess diagram: 

a{5r) = [0{g — r,r) — B{g — r, r)] ■ F{g — r, r + 5r) (2) 

and estimated the error of cross-correlation function as: 

a{a) = { [0{g - r, r) + B{g - r, r)] ■ F\g - r, r + 5r)Y-\ (3) 

g—r,r 

In the above equations, O and B represent the Hess diagrams of orbit and background 
segments, respectively. F is the value of M92 filter. 

Then we define the maximum position, ^r^, of the cross-correlation function to be the 
actual distance modulus to the stream. Near the maximum, we can estimate the cross- 
correlation function by the Taylor expansion: 

a{5r) = a{5rm) + (^) \5r=5r,„ 5r + ^(^) \sr=5r,„ 5r'^ + ... (4) 

The second term in the right equals zero, since the first derivative is zero at maximum. We 
use the above Taylor expansion, using only the lowest order non-zero terms, to estimate the 
error in the cross-correlation function: 

o-2(a) =< a^(5r) - a^(5r^) >= a('5'^m)(^) |5r=<5r,„< >^ . (5) 
Since < 6r >= asr, 



cr2(a) 



'^Sr — \l~r-. , J2„ \sr=5r^ (6) 

a(5r^)^ 
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The distance and corresponding error of each of the seven points along the stream 
with SEGUE spectra are hsted in Table 1. Figure 8 shows the same data as Figure 7, but 
overplotted with the M92 fiducial sequence shifted to the best estimated distance for each 
sky position. 



4. Notes on Individual Regions 

Figure 9 shows velocity histograms of low metallicity stars from Figure 6 within ±1.3° 
of each selected region. A Gaussian velocity histogram, representing a halo distribution with 
cr = 100 km s~^ and are normalized so that the area under the curve equals the number of 
stars in the histogram, is also shown. The region number is indicated in the upper right 
corner of each panel, and the velocity of the peak most likely associated with the GD-1 
stream is also indicated. We list below some details of each region's selection, including the 
SDSS/SEGUE plates on which most of the objects were detected. 

In each region with a clear stream detection, the velocity and velocity dispersion for 
the GD-1 stream were computed using an iterative method that used only stars within one 
standard deviation of the mean velocity. We computed the mean and standard deviation 
of the stars near the velocity peak. Then, we selected stars that were within one standard 
deviation of the mean and re-computed the mean and standard deviation. The standard 
deviation calculated this way is an underestimate, since we have removed the tails of the 
distribution. We corrected the standard deviation assuming Gaussian tails. This process was 
repeated until the computed mean and standard deviation matched the mean and standard 
deviation used to select the stars in the stream. 

Table 2 lists 48 high confidence GD-1 stream members. The sample includes the stars in 
Figure 6 (which are selected based on angular distance from the GD-1 stream, photometric 
color and magnitude, and metallicity) which have Galactic longitude within two degrees of 
the seven plate centers, have velocities near the measured or expected velocities for the GD- 
1 stream, and which have proper motions that are consistent with our GD-1 stream model 
(presented in §7). The velocity cuts for each of the seven regions are: 1) 97 < Vgsr < 118, 
2) 59 < Vgsr < 79, 3) 32 < Vgsr < 52, 4) -17 < Vgsr < 3, 5) -78 < Vgsr < -58, 6) 
—97 < Vgsr < —77, and 7) —166 < Vgsr < —146. The proper motions expected for each 
of the seven regions are listed in Table 1. The stars in Table 2 are within two sigma of 
the expected proper motions, where one sigma is 3 mas/yr in each of fii and fib- Table 2 
includes each objects SDSS-ID number, coordinates, magnitude, colors, velocity, estimated 
metallicity, surface gravity and proper n aotion. The prope r motions (/i/, /Xb) listed here 



are from the USNO-B/SDSS catalog of iMunn et all fl2004f ). as extracted from the DR7 
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database's 'propermotions' table. Errors on an individual (7 ~ 18 star's measurement are 
about 3 mas yr~^ on each coordinate. 



4.1. Primary Kinematic Regions 

Regions 4, 5, and 6 were specifically targeted by SEGUE with plates along the GD-1 
stream. These regions, along with region 1, targeted by the SDSS Legacy survey, constitute 
the 4 regions along the stream with spectroscopy used to fit a model orbit for GD-1. 

Region 4, Plates 2889 and 2914-' The GD-1 stars are well sampled in this plate pair; 
there is a clearly detected peak of (primarily) F turnoff stars at Vgsr = —7 ± 1 km s^^ with 
a = 3.9 km s^^, in the magnitude range 17.5 < go < 19.5. The metallicity of stars in the 
velocity peak is about [Fe/H] ~ —2.2. There are two BHBs in this sample near g ~ 14.9, 



corresponding to a distance of 7.2 kpc ( iSirko et al.ll2004l ). The distance to this region derived 



in Section 3 from the turnoff photometry is 7.5 ± 0.33kpc, in good agreement with the BHB 
magnitudes. The velocities of the two BHBs are -6 and -3 km s~^, in excellent agreement 
with the average of the turnoff star velocities, as are the metallicities ([Fe/H] = -2.0 and 
-2.1, respectively). 

Region 5, Plates 2557 and 2567: The plate spans the full width of the GD-1 stream. 
There is a strong peak in the velocity distribution. This peak has Vgsr = —71 ±2 km s~^ with 
a = 5.3 km s~^. The average g magnitude in this range is (7 ~ 19, which is consistent with 
the distance estimation from color magnitude turnoff fitting. The metallicity distribution of 
F turn-off stars in the peak shows that [Fe/H] of this stream is about -2.05. The distance to 
this piece of the stream is 8.0 ± 0.53 kpc. 

Region 6, Plates 2390 and 24IO: A strong and narrow peak is detected in Figure 9, 
at Vgsr = —87 ± 2 km s^^, a = 9 km s^^ and go = 18.82. [Fe/H] distributions show the 
metallicity peaks at about -2.05, which is consistent with the [Fe/H] peak find in plate 2567 
and slightly higher than plate 2914, this indicates that they are from the same stream. A 
distance estimate puts stars in this region of the stream at 8.8 ± 0.75 kpc from the sun. 

Region 1, Plate 1154: This is a special SDSS legacy plate, rather than a SEGUE plate, 
but since it is at low Galactic latitude it does have a large number of stellar candidates. 
Its important GD-1 stream candidate stars cluster at velocity Vgsr = 108 ± 5 km s~^, with 
(T„ = 11 km s~^, and (lower S/N) metallicity [Fe/H] = —2 ± 0.3, which anchors the stream 
orbit away from regions 4,5 and 6. With larger errors, the distance to this stream here is 
10.4 ± 1.2 kpc. There is a second peak (3 stars) at fgsr = —40 km s~^. We discount this 
second peak as being associated with GD-1, since two of its three members have [Fe/H] = 
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-1.8, significantly higher than the average for other stream peak members. This secondary 
peak could be related to the Anti- Center stream noted by iGrillmair. Carlin and Majewski 
(120081 ). as the radial velocity of these stars is RV ~ +112 km s~^, is close to their value for 
ACS-C (see lower panel of Figure 1 and Figure 2 of that work). 



4.2. Other Regions 

The following 3 regions were not involved in the model fit (see below), but as there are 
SEGUE spectra available here along the orbit defined by the imaging arc defined in Equation 
1, we examine these SEGUE plates for possible stream members. 

Region 2, Plates 1760, 2433, 2667 and 2671: This plate pair was targeted by SEGUE 
as it contains the open solar-metallicity cluster M67. The Sagittarius stream also passes 
near by, along with the Anti-center stream. By chance, the GD-1 stream arc appears to pass 
within 1.3° of the plate center, and several very low metallicity turnoff stars are detected 
with average Vgsr = 69 km s~^. The distance to the stream is 6.5 ± 0.64 kpc. 

Region 3, Plates 2304 (^^d 2319: Several stars with metallicity have velocities in a 
broad excess near Vgsr = 42 km s~^, however there is not a significant candidate peak here. 
A peak at fgsr ~ —105 km s~^ has stars with < [Fe/H] >~ —1.8, and is not a viable stream 
maximum. The distance from section 3 is estimated at 7.0 ± 0.36 kpc. 

Region 7, Plate 2539 and 2547: This plate pair is on the fitted arc of Equation 1, but the 
stream becomes too tenuous to be defined. We do not see a velocity peak here, but following 
the trend in the peaks of points 1-6 there may be two or three stars (above background) at 
'^gsr = —156 ± 10 km s~^, with the correct metallicity and proper motion to be associated 
with GD-1. We do not place as strong a confidence in the GD-1 membership of stars in 
this peak as the other 6 regions. A very uncertain distance estimate to the stream here, is 
about d = 9.9 ± 1.2 kpc from the Sun. Orbit fits, which were not fit to data in region 7, 
confirm that this is the correct velocity to be looking for stream stars, but that the distance 
is somewhat underestimated. 

We note that the estimated distances to individ ual regions are in very good agreement 



with those quoted by IGrillmair and Dionatod (120061 ) . indicating that the results are some 



what robust to details of the method and potential parameters. 
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5. The Observed Stream Properties 



The observed properties of the stars in the seven GD-1 stream candidate regions are 
summarized in Table 1, where we hst region number (N); Equatorial coordinates, {a, 6); the 
corresponding Galactic coordinates (/, b) with errors (we use 6 to denote a measured error, to 
distinguish it from the intrinsic dispersion, which we denote with the symbol a); the average 
Galactocentric standard of rest velocity and heliocentric radial velocity with an error, and 
the velocity dispersion. The velocity mean and dispersion were calculated as described in 
§3. The tabulated intrinsic dispersions are upper limits to the actual velocity dispersion of 
the stream; since they are similar in size to the velocity errors for each individual spectrum, 
the measurement is consistent with an intrinsic velocity distribution of zero. The error in 
the mean is the velocity dispersion divided by the square root of the number of stars used to 
compute it. Regions 2, 3, and 7 do not have clear, narrow peaks in the velocity distributions 
and therefore were not used to fit the orbit, though Figure 6 shows there are excess stars at 
about the right velocities. 

Table 1 lists the computed Galactic X, Y, Z positions (with respect to a right handed 
coordinate system with the Sun at (-8,0,0) kpc) for each stream region and a distance to 
the stream at that region, again with errors. For the four regions numbered 1,4,5 and 6, the 
velocity and position accuracies are highly significant, and there is high confidence of the 
GD-1 stream membership of stars highlighted in the corresponding boxes of Figure 6. 

Figure 10 highlights this confidence by showing a color-magnitude diagram of all high- 
confidence spectroscopic GD-1 stream candidates in regions 1-7, shifted to a standard dis- 
tance of 9 kpc based on the photometric Hess diagram fitting. These objects are also proper 
motion selected, in that only objects within ±2a of the proper motion of the best fit model 



are ke pt. Superimposed over the spectral candidates is the M92 fiducial sequence of lAn et al. 



(120081 ) shifted to a distance modulus of m — M = 14.76, identical to that used in Figure 1. 



Distances to the other regions were estimated as described in section 3. 

We now calculate our best metallicity estimate for GD-1 by selecting only the 48 stars 
with spectra in Figure 10. We note that these stars were pre-cut on metallicity at an 
earlier stage (Figure 6) to have —2.3 < [Fe/H] < —1.65. A histogram with bins similar to 
the measurement error yields a GD-1 stream metallicity of [Fe/H]= —2.1 ± 0.1 dex with a 
dispersion of a = 0.3 dex (essentially the measurement error). In addition to these statistical 
errors, the re may be systematic error s in the metallicity determinations from SDSS DR7 of 



0.2 dex flAUende Prieto et al.ll2008[ l. 
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6. Orbit Fitting 



We now use th e data listed in Table 1 to fit an orbit for the stream, assuming a fixed 
Galactic potential. iGrillmair and Dionatod (120061 ) postulated the progenitor of this stream 
is a globular cluster because it has a narrow width in the sky. As a cluster orbits the Galaxy, 
stars farther from the progenitor will depart from the orbit due to dynamical friction and 
scattering of the stream stars. Because the progenitor is presumed to be a compact object 
with a few km s~^ velocity dispersion, it is reasonable to assu me that the stars in th e 
tidal stream lie approximately on the orbit of the globular cluster (jOdenkirchen et al.ll2003l ). 
Dwarf galaxies, on the other hand, will experience larger spatial dispersions because they 
have larger dispersions in their energies. Therefore, in this paper, we fit the orbit to the 
positions and velocities of the stars in the tidal stream. 



6.1. Galactic Model 



The Galactic potential model used in this work follows directly from lLaw et. al.l (l2005l ). 
We use a 3 component potential - disk, bulge, and halo, of the form given in Equations ([7]), 
(ED, and (E]). 

GMdisk /„x 
^disk = -Q I = (7) 

'rI + (a + vz^ + b^y 

<^bul,e = (8) 

r + c 

^haio = vl^i. In (^Rl + ^ + (9) 

In these potentials, R^ = X'^ + Y^ and = X'^ + Y^ + Z'^, where (X, Y, Z) are Galactocentric 
Cartesian coordinates. We adopt the Sun - Galactic center distance so that Xsun = —8.0 kpc. 
The symbol Rc is the cylindrical radius from the center of the Galaxy, whereas R refers to 
the Sun-centered distance to an arbitrary point along the stream. In these potentials, M^isk 
and Mhuige are the masses of the disk and bulge, respectively. The spatial extent of the 
potentials are determined by a, the disk scale length, b the disk scale height, c the bulge 
scale radius, and d the dark matter halo scale length. Additionally, Vhaio describes the dark 
matter halo dispersion speed (which is related to the total dark matter halo mass), and q 
represents the dark matter halo fiattening in the Z direction. We found that the parameters 
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of the Galactic potential were not well constrained by the path of the ti dal stream, so the se 
parameters (Table [3]) were held constant, with the same values used by iLaw et. al.l (120051 ). 



In this work, we will fit four orbital parameters, given in Table HI (-R5, fa;,5, 'J^y.s, ^^2,5) 
are the distance (from the Sun) to, and velocities of, region 5. Given a Galactic potential, 
these parameters fully specify the GC orbit. We then evolve the test particle forward and 
backward from the s tarting locati on, using the mkorbit and orhint tools in the NEMO Stellar 
Dynamics Toolbox ( jTeubenlll995l ). We convert the resulting orbit into (Z, h) and perform the 
goodness of fit calculation. 

To find reasonable initial values for these parameters, we imagine placing a test particle 
in region 5 at (/,6, -R5) = (172.3°, 57.43°, 8.0 kpc). We then construct a vector between 
the (/,6,i?5) = (172.3°, 57.43°, 8.0 kpc) and (/,6,i?6) = (161.95°, 59.12°, 8.8 kpc) points. 
This gives us the direction of the total velocity, because we are assuming the orbit passes 
through both of these points. The principle initial values for the parameters in region 5 
are (-R5, fx,5, fy,5, ^2,5) = (8.0 kpc, -94 km/s, -285 km/s, -104 km/s). In practice we start 
searching for the best parameters in a range of values near these approximate values for the 
orbital parameters. The parameter selection ranges are given in Table HI 



6.2. Goodness-of-Fit 

In order to find the best-fit orbit to the data given in Table 1, it is necessary to define 
a goodness-of-fit metric, and search the relevant parameters for the minimum value of this 
metric. The metric for an orbit of this type involves three important factors: the orbit passing 
through the plate locations in the sky, having the proper velocities at these locations, and 
having the correct distances. In order to consider all of these factors, we define three chi- 
squared values, one for each of the relevant variables, and simply sum them to create the 
total goodness of fit metric. Specifically, 

2 \ ^ f ^ model, i bdata,i\ /-in,\ 

x, = l^[ J , (10) 

2 / '^'^ model A f'^dataA \ /-i -i \ 

^™ = z^ ' (11) 

Xl = (^ ^'"""'''^"^'"'"" ^ ' , and (12) 
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= - {xl + Xl + Xl) , (13) 

where t] = N — n — 1, N being the number of data points, and n being the number of 
parameters. 

To calculate these values, we calculate a model orbit using the selected parameters. 
We search the orbit for the I values from the plates, and use the associated values of b, Vgsr, 
and R to compute x^- 



6.3. Gradient Descent 

We now optimize the orbital parameters so that is minimized. To do this, we choose 
an initia l set o f parameters, calculate an orbit using the NEMO Stellar Dynamics Toolbox 
(ITeubenl Il995l ). and calculate x^- We then use a gradient descent method to adjust the 
parameters to new values, generate a new orbit, and recalculate x^- This process is continued 
until the true minimum value of is found, and the associated parameters describe the best 
fit orbit. 

Let the vector Q = (-R5, fz,5, ^^?;,5, ^^2,5) describe the parameters. For each Q there is an 
associated x'^iQ)- We choose an initial set of parameters Qo, and find x'^iQo)- We then 
iterate the parameters by Equation (ITlj) . 



Qi,new Qi,old hil\^ iX i,Q 



old 



(14) 



We calculate the gradient using a finite difference method, shown in Equation fllSI) . 



V. x\Q) 



X\Q^ + K)-X\Q^-hi) 



2h,, 



all other Qk fixed 



(15) 



Different values of hi are used because the parameters are on different scales, it would 
not be appropriate to use the same step size for them all. The step sizes for the parameters 
are given in Table HI 

A is a variable- learning parameter. It initially begins at A = 1, and if the new value of 
x'^iQ) is smaller than the old, then A is multiplied by 1.03, if not, it is multiplied by 0.80. 
The purpose of this is to ensure if a minimum is being found, then it is found faster than 
with a constant-learning parameter. We also multiply it by the associated hi value to make 
the step size appropriate for the parameter being considered. 
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6.4. Error Estimation 



The uncertainties of ttie best-fit parameters can be estimated from t he sfiape o f the 
surface at its minimum. To do this, we follow the method outlined by ICole et al.l (120081 ) . 
We construct a matrix V of second partial derivatives of the surface, evaluated at the 
minimum found by the gradient descent. The error estimate for the i^^ parameter is ai = 
\pWii. The V matrix is defined as 



V = H-\ 



(16) 



The matrix H is the Hessian Matrix, whose elements are given by 



H,, = '-^^-^ ^, where, (17) 

H]j = {Qj + hj ,Qi + hi)\ Q^j^g^ g^g^ 

H-j = x\Qj - hj, Qi + hi) I Q^j^g^ g^g^ 

= x\Qj + hj, Qi - hi) I Q^j^g^ g^g^ 

Htj = (<5j - /^i, <5.-/ii) I all other Qfc fixed 



7. Results and Discussion 

We select five initial sets of parameters by randomly choosing values within the ranges 
given in Table HI We perform the gradient descent to reach the best fit parameters. We 
then estimate the parameter errors using the Hessian method outlined above. The best-fit 
parameters and their errors are (-R5, fx,5, "f^y.s, "^^2,5) = (8.4 ± 0.8 kpc, —89 ± 2 km/s, —236 ± 
6 km/s, —115 ± 3 km/s). The chi-squared of this fit is 2.07. The negative velocities indicate 
a retrograde orbit. The perigalacticon for this orbit is located r = 14.43 ± 0.5 kpc from the 
Galactic center at (Z, b) = (158°, 60°), near region 6. The space velocity of stars in the model 
at this position is 276 km s~^. The apogalacticon is at r = 28.7 ± 2 kpc from the Galactic 
center, toward (/, 6) = (306°, —35°), though we do not observe this direction on the sky. All 
errors are la. 

These orbital parameters are fairly insensitive to the choices of parameters in the Galac- 
tic potential. In particular, they are good estimates of the orbital parameters for a wide range 
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of q and d. To investigate whether we could determine anything about the shape of the dark 
matter halo from this tidal tail, we performed parameter sweeps on q and d while keeping the 
kinematic parameters constant (Figure 11). We see that very low flattenings are inconsistent 
with the data, but a wide range of q and d are allowed. 

The potential assumed in fitting the orbit is a standard logarithmic flat-rotation curve 
dark matter halo plus a stellar disk. Since the GD-1 stream approaches within 15 kpc of the 
Galactic center, the effects of the massive disk are felt by the orbit, and increasing the relative 
mass of the disk vs. the halo can mimic the effects of a flattened halo. At perigalacticon, the 
disk exerts twice as much gravitational force as the halo. Given this, it is not surprising that 
the stream orbit depends minimally on the halo parameters. More models and constraints, 
from this and other streams, are clearly needed to constrain the shape of the dark matter 
halo. 

The upper panel of Figure 12 shows the orbit in (/, b) with the stream locations shown. 
The model prediction is in very good agreement with the experimental observations. The 
middle and lower panels of Figure 11 show the orbit in / versus fgsr and I versus distance 
from the Sun. We also see good agreement with the experimental observation. Figure 13 
shows the orbit projected into the three planes of Galactic coordinates (X, Y, Z). We deduce 
an orbital eccentricity e = 0.33 ± 0.02 (one sigma error) and an inclination to the Galactic 
plane of z ~ 35 ± 5°. Arrows show the relative direction of the stream's retrograde motion 
compared to the Milky Way. 

The final columns of Table 1 show the predicted proper motions (yU^, ^}^) for stars in 
the stream at each region 1-7 based on the distances in Table 1. These predictions may 
be compared with actual observed proper motions for stream star candidates in Table 2 at 
each region. In general the agreement is quite good for regions 2-6, given proper motion 
errors of Icr = 3 mas yr~^ in each coordinate. Spectral candidates more than 2a away were 
excluded from Figure 10, dropping about 20% of the candidates, leaving a generally good fit 
to a shifted M92 fiducial sequence for these regions. Region 7 had fewer good proper motion 
matches, and it is possible that we are not seeing GD-1 stream candidates here. 

Figure 14 shows the photometrically chosen stars with proper motions available near 
regions 1,4,5 and 6, along with an equivalent set of field stars (chosen 5 degrees away) for 
comparison. There's a clear excess of 'on-stream' stars in the lower right quadrant of each 
region on-stream - these are likely stream members. The estimated tangential velocities 
(relative to the Sun) are given. 

To search for a possible progenitor, we selected all Milky Way globular clusters from 



Harris! (119961 ) that had metallicities in the range —2.5 <[Fe/H]< —1.5. Only seven of these 
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globular clusters (Terzan 8, Arp2, NGC 6809, NGC 6749, NGC 6341, NGC 6681, and NGC 
6752) are within 5° of the GD-1 orbit. Additionally, we considered NGC 2298, which is 5.5° 
from the orbit, and has a metallicity of [Fe/H] = —1.85. We compared the positions and 
velocities of these globular clusters with an orbital path that extends all the way around the 
Milky Way. To create a stream of this length would require a globular cluster to orbit the 
Milky Way for on the order of gigayears, with the length depending on the concentration 
of the progenitor, as well as the shape and location of the progenitor's orbit. Of the eight 
globular clusters, NGC 6809, NGC 6749, and NGC 6752 are ruled out because their distances 
are more than a factor of two different from the distance to the orbit. The remaining five 
clusters had radial velocities that are inconsistent with the predicted orbit by more than 
50 km/s. We therefore conclude that the Milky Way globular cluster catalog published by 



Harris! (119961 ) does not contain the progenitor of this stream. 



8. Conclusions 

We use spectroscopic kinematic and abundance information to isolate stars in the GD-1 
stream, and use the positions and velocities of those stars to derive orbital parameters for 
its orbit. The GD-1 stream is moving very rapidly on a retrograde orbit around the Milky 
Way. In the region of the orbit which is detected, it has a distance of about 7-11 kpc from 
the Sun. The stream's orbit takes it to apogalactic distances of 28.75 ± 2 kpc, and it has a 
perigalacticon of 14.43 ± 0.5 kpc, implying an eccentricity of 0.33 ± 0.02. The inclination to 
the Galactic plane is about i ~ 35° ± 5. The metallicity of the stream is [Fe/H]~ —2.1 ± 0.1 
plus systematic errors of a few tenths dex. None of the known globular clusters in the Milky 
Way have positions, radial velocities, and metallicities that are consistent with being the 
progenitor of the GD-1 stream. 

The consistency between the proper motions of these stream candidates and our best 
fit model gives us further confidence that we have identified stream members and that our 
model accurately represents the path on the sky of the stream stars. While we claim only 
consistency here between the proper motion data and our model, we note that more detailed 
fits to the proper motion (in addition to the radial velocities) for such nearby streams can 
be a crucial tool in constraining the halo potential shape and other parameters. 
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Fig. 1. — The GD-1 stream Hess diagram 

Hess diagram of on-stream minus off-stream (taken from 5 degrees lower 6 than the 
on-stream stars) star counts. Darker areas of the figure have higher stellar density. Thick 
disk stars and M dwarfs did not subtract perfectly between on and off regions. The blue 
turnoff at {{g — T)o,fi'o) = (0.25, 19) indicates that t he GD-1 stream is relatively met a l poor . 



A fiducial sequence for globular cluster M92 from IClem. Vanden Berg, fc StetsonI (l2008l ) , 
shifted to m — M = 14.76, is overlaid with black dots. The region of color-magnitude space 
used to select the photometry of stars in Figure 2 is outlined in heavy black lines. The region 
of color-magnitude space that is used to select the stars plotted in Figure 3 is outlined in 
light black lines. 

Fig. 2. — The GD-1 stream photometry 

The upper panel shows a Hess diagram of stars in the SDSS footprint which are within 
the color and magnitude boxes thickly outlined in Figure 1. The GD-1 stream arcs faintly 
from [a, 6) = (220°, 58°) to (126°, 0°). Seven regions where spectroscopy of GD-1 stream star 
candidates have been obtained are numbered. Squares indicate regions where the stream was 
clearly found in velocity, and circles indicate additional plates that may contain stream stars. 
The lower panel shows the GD-1 stream in Galactic polar coordinates. The color-magnitude 
cuts are similar to the above, but plotted in Galactic polar coordinates as labeled. Note the 
stream's path on the sky. 

Fig. 3. — Velocities of stars near the stream 

We plot the line-of-sight, Galactocentric standard of rest velocity versus Galactic lon- 
gitude for SEGUE and SDSS stars for which we have spectra, and which are close in color- 
magnitude space to the fiducial sequence of M92 (as indicated by the narrow dark outlines 
in Figure 1), which are close (within 1.3°) in projected distance to the GD-1 stream shown 
in Figure 2. A sine curve with amplitude 110 km s~^ is shown to indicate the locus of 
stars rotating with the Sun about the Galactic center. Stars in the halo will have velocities 
centered on Vgsr = and a large a ~ 100 km s~^ dispersion. Regions where GD-1 stream 
candidates are followed up on are numbered 1-7. Note in particular the groups of stars at 



V. 



gsr 



-90 km s ^ in regions 5 and 6. 



Fig. 4. — Histogram of Vgsr towards regions 5, 6 

All stars with spectra from Figure 3 towards regions 5, 6 of Figure 2 are histogrammed in 
Galactocentric velocity. Gaussians are overlaid representing a thick disk and halo distribution 
toward this direction. The candidate GD-1 stream stars at about Vgsr ~ "82 km s~^ cannot 
be explained by either a halo or thick disk distribution. 
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Fig. 5. — Metallicity distribution of region 5,6 stars 

Histogram of SSPP metallicities for all stars (light line) and for stars with —97 < Vgsr < 
—61 km s~^ (heavy hne). The velocity-selected peak has an excess of stars with metaUicity 
lower than that of the halo. The hashed line indicates a correction to the heavy line for 
interlopers at other velocities which have the same metallicities as candidate stream stars. 
This figure suggests that the peak metaUicity of the GD-1 stars is lower than [Fe/H]=-1.9. 

Fig. 6. — Velocity distribution of low metallicity stars near the stream 

The subset of stars from Figure 3 with SSPP metallicities —2.3 < [Fe/H] < —1.65 are 
presented. Now the counter-rotating GD-1 stream stars stand out more clearly against their 
field contaminants compared with Figure 3. Regions 1,4,5 and 6 have clear peaks; the mean 
and the error on the mean for these peaks are indicated by the position and height of the 
rectangles at these four longitude locations. The circles in regions 2, 3, and 7 indicate the 
area through which the stream should pass if our model is correct. All areas except region 
7 seem to have an excess of stars with the expected stream velocities. 

Fig. 7. — Background Subtracted Hess diagrams along the stream 

At each region 1-7, stars with {u — g)o > 0.4 within 0.5° of the cubic described by 
equation 1 are selected. A similar background set of stars, obtained from a region 1.5° away 
from the cubic, was selected for each of the seven regions. The Hess diagrams show the 
stream star counts minus the background star counts, as a function of go and {g — r)o. Light 
(white) regions indicate an excess of stars along the stream. 

Fig. 8. — Hess diagrams with fits 

The diagrams from the previous figure are fitted with M92 fiducial loci, shown as (green) 
crosses, shifted up and down in magnitude until a maximum correlation is obtained (see text). 
Candidate stream stars with spectra are indicated as (red) diamonds. The distance to each 
region along the stream is then derived from the best fitting shifted M92 fiducial. 

Fig. 9. — Velocity selection of stream candidates 

Velocity histograms in each GD-1 stream candidate regions 1-7, with SSPP metallicities 
—2.3 < [Fe/H] < —1.65. A Gaussian with a — 100 km s~^ is overlaid. The Vgsr of the 
expected GD-1 peak is noted in each panel. 

Fig. 10. — Stars with spectra in the stream 

All stars in regions 1-7 with SEGUE spectra which meet the metallicity, color, mag- 
nitude, velocity and proper motion cuts as described in the text are plotted with colored 
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points as indicated in the legend, along with estimated distance to each set of points. Each 
set of points was shifted to the reference distance of 9 kpc (of Figure 1), and overlaid with 
a M92 fiducial locus, shifted to the same distance moduli. Note the two (blue) points at 
5^0 ~ 15.3, which are actually at go ~ 14.9 before shifting. These are candidate BHB stars 
in the GD-1 stream, at a distance of 7.5 kpc from the Sun in region 4. 

Fig. 11. — Parameter Sweep of q and d 

Fixing the orbital parameters at their best fit values, derived assuming that q — 1 and 
d = 12, we varied d and q and measured x^. This is justified by our observation that the 
best fit GD-1 orbital parameters are fairly similar for a large range of assumed q values. The 
normalization of this surface in d, q was chosen so that the {d, q) — (12, 1) is 2, to 

match the optimization presented in this paper. The dark squares in the figure show the 
locus of lowest in this two dimensional surface. We see that very low values of q are not 
favored, while there is almost no constraint on how prolate the dark matter distribution can 
be. The values of q and d are coupled. The best fit d, for a given q shifts somewhat with 
small changes to the orbital parameters, so this plot cannot be used to definitively pin down 
the d value, even though the x^ surface seems narrow in that dimension. 

Fig. 12.— Best Fit Orbit 

Upper panel: Galactic coordinates. The best fit orbit (with fixed g, d) to the data in the 
four regions 1,4,5 and 6. Region 1 is leftmost, with I = 224.47°. Galactic (/,6) for all seven 
regions described in the text are plotted as crosses. Middle panel, plotting Galactic longitude 
/ vs. Galactocentric radial velocity Vgsr for the best fit model and data. Regions 1,4,5 and 
6 (left-to-right) have the smallest error bars. Lower panel, plotting Galactic longitude I vs. 
Sun-centered distance to the stream. The errors on regions 4, 5 and 6 are small; the error 
bars on the other region data points are limited by how well the position of the stream turnoff 
(minus a background field) can be identified in a color-magnitude Hess density diagram of 
stream stars. The distance error bars are largest for region 7 (rightmost), where the stream 
is confusion limited. 

Fig. 13.— Best Fit Orbit: X, Y, Z 

Shown in heavy black line is the best fit orbit to the four regions 1,4,5 and 6 in Galactic 
rectangular coordinates {X, Y, Z) in each of the cardinal projections. The coordinate system 
is a right-handed system with the Sun at at (-8,0,0) kpc and the Galactic center at the 
origin. The coordinate and turnoff magnitude data from the seven regions is converted to 
{X, Y, Z) assuming an absolute turnoff F star magnitude of Mg — -1-4.2. The seven points 
are indicated with error bars from distance error estimates. Regions 1 and 7 are indicated 
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in each panel, with the other regions faUing in order at intermediate positions. The units 
of each axis is kpc. An arrow originating at the Sun indicates the direction of Galactic 
rotation. The arrows associated with the stream indicate the retrograde direction of motion 
of GD-1 stream stars. The space velocity of the stream at perigalacticon is approximately 
276 km s-^ 

Fig. 14. — Proper Motion consistency 

For the four regions N=l,4,5,6 where we have fitted the orbit, we select stars within 
0.3° of a cubic similar to equation 1 which have measured USNO-B proper motions from 
the DR7 database. For reference, we select similar sets of field objects offset by 5 degrees in 
declination from the on-stream objects. We then sub-select stars with colors and magnitudes 
of stream turnoff candidates using the the heavily outlined box of Figure 1, and plot the /x; 
vs fib of the on-stream (black dots) and off-stream (open circles) for each selected region. 
There is a clear excess of on-stream points extending to the lower right in each region. We 
superimpose crosses representing the point where the best fit model (Table 1) crosses the 
stream, and show the derived tangential velocity in the Figure. The cross always falls in the 
same quadrant with the excess of proper motion points. Typical errors on each point are 
3 mas yr~^ in each direction. The upper and right axes in each figure convert the observed 
proper motions to tangential velocities, assuming the distance to the stars is the distance to 
the fitted orbit for that particular GD-1 stream region. 
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Table 1. GD-1 Stream Detections 
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Table 2. Spectra of GD-1 Stream Candidates 



SpecID 


N 


a 


5 


l?gsr 


a(v) 


RV 




{u-9)o 


{g - r)o 


[Fc/H] 


^([Fc/H]) 


log g 


cr{logg) 
















km 


km 


km 


mag 


mag 


mag 


dex 


dex 


dex 


dex 


mas / yr 


mas/yr 


li04-Ov5Uoo-ZOD 


1 


IzO.Doooy 1 


-u.4oyooi 


1 1 Q /I 
1 1 • 4 


7 A 
i .4 


Ofi/1 7 
Z04. / 


1 O 1 07 

ly . lu ( 


O OOO 

u.yuy 


n 97f^ 
U.z i 


1 ftO 

-i.oy 


O Oft 

U.Uo 


Q Afi 
O.40 


o 9c; 
U.zO 


Q 

O.O 


A Q 
-4.0 


li04-0oUoo-i40 


1 


IzD.O 1 ( lOU 


O /1QO'3Qfi 

-u.4oyooo 


OV 

y / .o 


7 Q 


0/1Q /I 
Z4y.4 


1 Q 1 oo 

ly. luz 


O 00/1 

u.yu4 


n 1 fto 
U. loU 


9 1/1 
-Z. 14 


O 07 
U.U ( 


Q *^0 


O 9Q 

U.ZO 


1 9 Q 

iz.y 


1 O 
-lU.O 


ll04-0oUO0-l00 


1 
1 


1ZD.D404U ( 


O 'i.ACxATI 
-U.o4U4 ( / 


1 OQ O 

luy.u 


Q 

O.O 


O^IO Q 
zOU.o 


1 ft 1 1 7 
lo.ll ( 


O QQC; 

u.yoo 


o 97c: 
U.z (0 


9 O^l 
-z.UO 


o riA 

U.U4 


Q QO 

o.oU 


O 99 
U.ZZ 


1 1 ^; 

11.0 


Q Q 

-O.y 


i /DU-OoUoD-ooy 


z 


lol.4UDycS4 


O '7'70/1 1 Q 
y. I ( U4l0 


DO.fS 


1 /I K 
14.0 


loo.y 


1 ft Qftn 

lo.ooy 


O OOO 

u.yuy 


O 007 

U.zz { 


O 1 ft 

-z.lo 


oil 


A OK 

4.U0 


o oo 

u.zy 


O A 

y.4 


ft A 

-0.4 


ZD 1 1-04 14 i- 40 Z 


z 


1 Q9 1 

loz.io 1 cSOo 


111 c^Qc;i 
ll.loOttOl 


/ (.o 


O.U 


1 QQ Q 

lyo.y 


1 ft 7ftO 

lo. / oU 


O ft/11 
U.o41 


n Ana 
U.4U0 


O 0/1 
-Z.Z4 


O HA 
U.U4 


Q KA 
O.04 


O 90 
U.ZZ 


1 o o 
lU.U 


c; o 
-O.z 


ZOD (-0414Z-4Z / 


z 


loz.oO f 000 


I 1 Q9C^/1QC; 

II .ozo4yo 


70 c; 
( U.O 


A n 

4.U 


loD.O 


1 7 y1 QQ 
1 ( .4oo 


fl QS7 

u.yo ( 


o /loo 
U.4zU 


9 17 
-Z.l / 


O OO 
U.Uz 


9 77 
Z. i ( 


O 9/1 
U.Z4 


(\ 7 
0. / 


7 A 
- / .4 


ZOD /-0414Z-DU4 


z 


loo.U 1 ooyz 


IZ. ZD0oo4 


f;9 Q 

DZ.O 


9 K 
Z.O 


1 7c; 1 
1/0.1 


lo.ooy 


1 1 f^O 

i.ioy 


O A ftfi 
U.4oO 


O Oft 

-z.Uo 


O OQ 

U.Uo 


9 99 
Z.ZZ 


O 1 
U. 10 


ft A 
0.4 


d 1 
-0.1 


0Q1O f^Q^KQ Q/IV 


Q 

o 


loo.oyyo ( 


99 A 1 Q099 


QQ 1 

oo. 1 


1 1 fl 
11. D 


114.0 


1 O /I f^l 

ly .401 


O ft/1 Q 
U.o4o 


n 9/1 Q 

U.z4v5 


O 0/1 
-Z. Z4 


o oo 

U.ZZ 


Q ftO 


O OO 

u.uy 


1 Q 9 
lO.Z 


7 7 


Zoiy-OO IDO-ODO 


Q 

o 


1OO.40U ( 01 


99 /I Q/1 f^79 
ZZ.4O40 1 Z 


/10 /I 
4U.4 


/f 7 


1 1 R 7 
110. / 


1 7 OOfi 

1 / .yyo 


U.oOO 


n 9c;o 

u.zoy 


1 OQ 


O OO 
U.Uz 


Q OO 


O 1 o 
U. lU 


1 9 
IZ.O 


O 7 

-y. i 


ZOiy-OO ( DO-OO / 


Q 

o 


lOO.DOyo ( O 


zz. ooy / oo 


/I Q O 

4o.y 


1 1 Q 
11. 


1 0/1 7 

lz4. / 


1 O fiOQ 

ly.ouo 


O 07/1 

U.O / 4 


o Qftc; 
U.ooO 


O 07 

-z.U / 


O 07 

U.U ( 


Q dQ 

o.Oo 


O QO 

U.oU 


1 O 1 

Iz.l 


7 O 

- / .y 


ZOc5y-040oU-oll 


4 


1 AO VfiVQI c; 
14Z. ( O ( olO 


OQ A«^ nac^ 
zy.401UoO 


n Q 
-u.y 


A O 
4.Z 


/t7 ft 
4/ .o 


1 ( .404 


O Q70 

U.y ( z 


0/117 
U.41 ( 


1 Oft 

-i.yo 


O 0/1 
U.U4 


Q OQ 

o.zo 


O 9ft 

U.ZO 


1 O 7 
lU. ( 


Q A 

-y.4 


OOI /I C^/IC^QQ 00*7 

zyi4-o4ooo-zy / 


/I 
4 


14z.y / oc5U4 


oo 1 QfiQO'T 

zy.ioooU / 


Q O 

-o.U 


Q n 

o.y 


/I n 
40. y 


1 ft 1 QO 

lo. loy 


O OOC^ 

u.yzo 


O 1 ftO 

u.loy 


o oo 

-z.uy 


o r\A 
U.U4 


Q nA 

o.y4 


O 1 

U.IO 


ft ft 
O.O 


O Q 

-y.o 


zooy-040oU-z4U 


A 
4 


1 /IQ Q "iOC^O/l 
14o.OOyOZ4 


90 fl90'iOQ 

zy.Dzuoyo 


9 

-O.z 


A 7 
4. / 


/1 1 7 
41. / 


1 7 QQQ 

1 / .yoo 


O ftftQ 

U.OOO 


U.zOO 


9 1/1 
-Z. 14 


O OQ 

U.Uo 


Q KQ 

O.Oo 


O 1 7 
U.l / 


O 

y.o 


a 7 

-0. ( 


Zooy-040oU-ZU4 


/I 


1 /I Q /1 9/1 DOS 


Zy .UOD40O 


-O.O 


^ 1 
0.1 


4D.0 


1 7 899 
1 / .oZZ 


n ft79 

U.O i z 


n 971 

u.z i 1 


9 91 
-Z. Zl 


U.UO 


9ft 
o. Zo 


n 90 

U.ZU 


(\ 7 


A <^ 
-4.0 


9QCQ KAK'in 91 
Zf5oy-040oU-ZlO 


4 


^ A'i A c;qi oq 
14o.40olof5 


OQ 1 90fi'i9 

zy. izuDoz 


Q /I 
-0.4 


1 7 
1. ( 


AR A 
40.4 


1 A Qfti 
14.yol 


1 1 Q/1 

1. iy4 


O 9Q7 
-U.ZO ( 


9 OQ 

-Z.uy 


oil 

u.ll 


Q OO 

o.zU 


O 1 7 

u.l / 


Q 

y.o 


1 O Q 

-lu.y 


ZOoy-040oU-Z0O 


4 


14o.00zy4c5 


zy.4( fouo 


-4.0 


A K 
4.0 


AA 1 
44. 1 


1 7 H/1 Q 

1 r .o4o 


O ftftft 
U.OOO 


n OQQ 

u.zyo 


9 01 
-Z.Zl 


O 01 

u.ui 


Q an 
O.OU 


n 1 ft 

U.lo 


ft A 
0.4 


Q 1 

-y.i 


Zoc5y-040oU-ZZ0 


4 


14o.0y 1 ooO 


OQ QOOAQ? 

zy.ouzoy ( 


R 1 
-D.l 


1 7 
1. ( 


A^ C\ 
41. U 


1 A ft/1 
14.O40 


1 1 Qd 
i.iyo 


o 1 oo 
-U.lUz 


1 Qft 

-i.yo 


U.UO 


Q Q/1 
o.o4 


O OQ 

u.uy 


1 Q 7 
lo. ( 


111 
-11.1 


9Q1 /I C^QQ 1 71 

zyi4-040oo-i f 1 


A 
4 


1 /I Q Vl QQK9 

i4o. f loyoz 


OQ Q.TA 1 on 
zy.o / 41UU 


-O.O 


o.U 


/I Q O 
40. Z 


1 Q /I QO 

iy.4oz 


n 7/17 

U. f 4 / 


o Oftc; 
U.zoO 


9 OA 

-z.zo 


O Oft 

U.Uo 


Q QQ 
O.OO 


o 

U.OD 


^^ A 

10.4 


ft 7 
-o. 1 


zy i4-04000-40o 


4 


1/1/1 1 07*5/1 
144. lU / O40 


QO 7Q700/1 
oU. ( o ( yz4 


O Q 

-y.o 


/I s 
4.0 


QQ ^1 
OO.O 


1 ft 71 
lo. i 10 


O 770 
U. 1 (u 


O 91 fi 

U.ZIO 


O Ofi 

-Z.zo 


o oc; 

U.UO 


A Al 
4.41 


O 1 o 

U. lU 


1 Q 1 
lo. 1 


-0.0 


zyi4-o4ooo-4oo 


/I 
4 


1/1/1 1 OOO/1 1 
144. lzUy41 


QO 0/1 770/1 

v>U.y4 1 ( z4 


1 A 
-1.4 


A A 
4.4 


/II O 

41. z 


1 7 Ot; K 

1 / .yoo 


O ft07 

U.oz / 


O OOO 

U.zzU 


O 1 o 
-z.lU 


o oo 
U.Uz 


Q 70 
O. IV 


O Oft 

U.ZO 


1 /I Q 
14. 


fl Q 

-O.O 


901 A t^A C^'^Q l^nO 

zy 14-040 oo-ouy 


/I 
4 


1/1/1 9'7c;'790 
144. Z 1 / ZU 


QO OfiO/H^O 

oU.zoU4Do 


A A 
-4.4 


Q O 

o.y 


/lO 
4U.0 


1 ft OOO 


O OOO 

u.yuy 


O 1 ft7 

U. lo / 


O Of^ 

-Z.zo 


o oi 
U.UI 


Q /Ifi 


O Q1 

U.ol 


o 

y.o 


1 O O 
-IZ.Z 


zy 1^-04000-010 


/I 


1 AA Aii70'iA 
I'i^.'lO ( Z04 


QO /II O/IOI 
OU.41U4U1 


A 7 


O.O 


"XGk ft 

oy.o 


1 ft 9/17 


U.oOO 


u.zoz 


1 ft9 

-1 .OZ 


n c\A 

U.U4 


79 
O. ( Z 


U.IO 


1 n 

lU.O 


1 

-l.D 


2914-5453O-540 


4 


144.561728 


30.968385 


-1.9 


11.0 


40.4 


20.190 


0.854 


0.344 


-1.95 


0.10 


3.46 


0.21 


12.6 


-3.8 


2567-54179-396 


5 


157.181542 


44.514349 


-67.9 


6.8 


-85.0 


18.745 


0.841 


0.254 


-2.08 


0.04 


3.78 


0.19 


13.9 


1.8 


2567-54179-246 


5 


157.565788 


43.643316 


-73.0 


6.0 


-87.1 


18.630 


0.934 


0.167 


-2.23 


0.01 


3.71 


0.18 


12.9 


-0.9 


2567-54179-238 


5 


157.691771 


43.797748 


-71.9 


12.7 


-86.7 


19.546 


0.908 


0.206 


-2.21 


0.02 


3.27 


0.16 


10.1 


-5.6 


2567-54179-189 


5 


157.817344 


43.923322 


-70.1 


11.9 


-85.5 


19.761 


0.772 


0.342 


-1.88 


0.07 


4.23 


0.36 


9.0 


3.5 


2567-54179-211 


5 


157.889540 


43.500286 


-65.2 


10.0 


-79.0 


19.309 


0.901 


0.205 


-1.89 


0.15 


3.74 


0.07 


9.2 


-0.3 


2567-54179-212 


5 


158.051293 


43.644022 


-67.6 


5.2 


-82.1 


18.349 


0.839 


0.208 


-1.90 


0.07 


3.67 


0.26 


9.7 


-0.2 


2567-54179-170 


5 


158.137638 


44.086320 


-58.2 


9.8 


-74.4 


19.112 


0.863 


0.201 


-2.04 


0.06 


3.64 


0.17 


14.4 


-1.7 


2567-54179-491 


5 


158.437074 


44.466389 


-69.2 


9.1 


-87.1 


18.242 


0.880 


0.209 


-2.02 


0.09 


3.95 


0.16 


8.9 


-1.9 


2567-54179-511 


5 


158.700621 


44.530392 


-68.9 


4.4 


-87.2 


18.083 


0.954 


0.237 


-2.22 


0.01 


3.18 


0.25 


12.1 


1.3 


2557-54178-498 


5 


158.747249 


44.686776 


-74.2 


3.9 


-93.2 


16.794 


1.048 


0.476 


-2.12 


0.02 


2.60 


0.27 


13.5 


2.7 


2410-54087-297 


6 


162.371533 


47.207633 


-87.0 


6.9 


-118.3 


18.893 


0.885 


0.255 


-2.26 


0.02 


3.46 


0.20 


15.4 


-1.3 



Table 2 — Continued 
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Table 3. Fixed Galactic Potential Parameter Values 



Parameter 


Value 


a 


1.0 




1.0 X 1O"M0 




3.4 X IO^^Mq 


a 


6.5 kpc 


b 


0.26 kpc 


c 


0.7 kpc 


'"halo 


114 km/s 


g 


1.0 


d 


12 kpc 
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Table 4. Fit Parameters 



Fit Parameter 


Description 


Step Size 


Random Selection Range 




Sun-centered distance of Region 5 


0.1 kpc 


5 : 13 kpc 


Vx,b 


X Velocity of Region 5 


1 km/s 


-130 : -60 km/s 


Vy,5 


Y Velocity of Region 5 


1 km/s 


-320 : -230 km/s 


Vz,5 


Z Velocity of Region 5 


1 km/s 


-80 : -120 km/s 
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